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Continuum coupling correction to binding energies in the neutron rich oxygen and fluorine isotopes 
is studied using the Shell Model Embedded in the Continuum. We discuss the importance of different 
effects, such as the position of one-neutron emission threshold, the effective interaction or the number 
of valence particles on the magnitude of this correction. 
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Nuclear structure of light neutron-rich nuclei has at¬ 
tracted much attention because of unusual spatial fea¬ 
tures of these nuclei and the isospin effects in the 
magic numbers |^. In this respect, the systematics of 
neutron separation energies Sn and interaction cross- 
section in sd shell region seem to indicate the disappear¬ 
ance of = 20 closure and the appearance of new shell 
closure at N = 16 in neutron-rich nuclei with Tz > 3 
close to the drip-line ||]. Another astonishing property 
is a large change of the neutron drip-line by adding one 
proton to the oxygen core, as indicate the instability of 
oxygen isotopes for TV > 16 |^,|| and the stability of 
with N = 22 ||] . Theoretical description of the low-lying 
states properties in these nuclei requires taking into ac¬ 
count the coupling between discrete and continuum cor¬ 
related states . This aspect is particularly important 
near drip-line where one has to use both structure and 
reaction data to understand basic shell structure and dy¬ 
namics of these nuclei ||,^. In weakly bound light nu¬ 
clei, not only coupling of bound states to continuum is 
strong but also many-body correlations should be treated 
rigorously. Both the continuum coupling and the realis¬ 
tic configuration mixing can be taken into account in 
the Shell Model Embedded in the Continuum (SMEC) 
1^. This approach provides a unified description of di¬ 
verse nuclear characteristics such as the energy spectra, 
including nucleon emission widths and electromagnetic 
transition probabilities, and the reactions involving one- 
nucleon in the continuum . In this work, we analyze 
in the Shell Model (SM) framework the continuum cor¬ 
rection to binding energies and neutron separation ener¬ 
gies. 

The detailed description of SMEC formalism has been 
given elsewhere ^,0. Separation of (quasi-) bound {Q 
subspace) and scattering {P subspace) states is achieved 
using the projection operator technique Q. P contains 
asymptotic channels made of {A — l)-particle localized 
states and one nucleon in the scattering state. Local¬ 
ized many-body states which are build up by the bound 
and resonance single-particle (s.p.) wave functions, are 
included in Q. They are obtained by solving standard 
multiconfigurational SM problem for the Hamiltonian 


Hqq. The residual coupling Hpq between states in Q 
and P is given by [D : El 2 = [a + (3P^2 ]^{yi - r 2 ), 

where a j3 = \ and is the spin exchange opera¬ 
tor. The effective SM Hamiltonian including the cou¬ 
pling to the continuum is energy-dependent : H^q{E) = 

Hqq + HqpG^p\e)Hpq, where G^p\E) is a Green 
function for the motion of single nucleon in the P sub¬ 
space. H^q is a complex-symmetric matrix above 
and Hermitian below it. The energy scale is settled by 
the one-nucleon emission threshold if |P. Radial s.p. 
wave functions in Q and the scattering wave functions 
in P are generated in a self-consistent procedure, start¬ 
ing with the average potential of Woods-Saxon type with 
the spin-orbit and Coulomb parts included, and taking 
into account the residual coupling. This procedure yields 
new orthonormalized wave functions in Q and P and new 
self-consistent potentials for each many-body state in Q. 

In the present studies we use the full sd valence space 
for N < 20 and the full pf shell for N > 20. For the 
effective interaction in Hqq we take USD Hamiltonian 
for the sd shell and the KB interaction for the pf 
shell . The cross-shell interaction is the G-matrix . 
In the calculation of binding energies from the nuclear 
energies in SMEC and SM, we proceed as in Refs. |12|. 

The ground state (g.s.) continuum coupling correc¬ 
tion in SMEC is calculated as : Ecorr =< ‘J’g.s.l'ffgQ — 
HqqI^q.s. >, the g.s. wave function in parent nucleus 
{N, Z) is coupled to different channel wave functions, 
which are determined by the motion of an unbound neu¬ 
tron relative to the daughter nucleus {N —1,Z) in a cer¬ 
tain SM state The magnitude of coupling matrix 

elements vary depending on the structure of SM wave 
function and the value of one-neutron emission 


threshold En'^'^K The individual contributions for differ¬ 
ent channel couplings tend to decrease with increasing 
excitation energy of a SM state due to a growing 

mismatch of both quantum numbers and radial proper¬ 
ties of wave functions in parent and daughter nuclei . 
This decrease is exponential on average. In the present 
studies, we include all asymptotic channels composed of 
SM states; their number vary from 3 in to 1837 
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in Even though in most nuclei these channels are 
closed, nevertheless the coupling to them modifies the 
binding energy depending both on and the struc¬ 

ture of the g.s. wave function. 




FIG. 1. The systematics of Ohui SMEC correction in the 
oxygen isotopes, (a) Neutron number dependence of the 
SMEC energy correction to the SM g.s. energy. The solid line 
is obtained for one-neutron emission threshold calcn- 

lated in SMEC for each nucleus. The dotted line with squares 
and triangles is obtained for which is fixed arbitrarily 

at 0 and 4 MeV, respectively, (b) One-neutron separation en¬ 
ergies Sn calculated in SM (dashed line) and in SMEC (solid 
line) are compared with the data (closed circles with the er¬ 
ror bars) 0. (c) Two-neutron separation energies S 2 n in SM 
and SMEC. Symbols are the same as in (b). 

The choice of the initial average potential follows a 
procedure described in Ref. 0. We demand, that the 
self-consistent potentials yield the energies of s.p. states 
at the same position as in the average potential with 


Woods-Saxon central part of radius Rq = fm, 

depth Vo{N,Z) = T 0.69(fV - Z)/{N + Z)] 

MeV {—/+ stands for neutron/proton potential), dif¬ 
fuseness a = 0.67 fm. The spin-orbit potential is 
ysoA^(21 • s)r-^dfir)/dr with A" = 2fm2, V,oiN,Z) = 
0.131Vb(-A^) Z) MeV and radial formfactor /(r) of the cen¬ 
tral potential. The Coulomb potential is calculated for 
the uniformly charged sphere with radius Rq. 

Fig. 1 shows results for oxygen isotopes. In our 
Ohuj description, the continuum coupling in these nuclei 
contains only neutron-neutron (T = 1) part. For the 
strength — (3) = 414 MeV-fm^, one ob¬ 

tains a good overall agreement with the spectra of oxy¬ 
gen isotopes. For example, SMEC calculation in 
yields the first excited state 2^ above the experimen¬ 
tal En^’^\ in agreement with the data |0]. The nn- 
continuum coupling modifies T = 1 monopole terms |l^ 
of the SM interaction. A good agreement for the binding 
systematics in oxygen chain can be found provided the 
T = 1 part of two-body interaction is modified as follows 


: < 5143=1 


= -0.05 MeV, 


= - 0-35 MeV. 


1/2 i0<i3/2 


= -bO.25 MeV, 


Fig. la shows the neutron number dependence of 
Ecorr for (i) En^’^^ of SMEC (solid line), and for (ii) 
En^^'^ fixed arbitrarily at 0 and 4 MeV (the dotted lines 
with open squares and triangles, respectively). The N- 
dependence of Ecorr in these three cases exhibits approx¬ 
imately quadratic dependence on the number of valence 
neutrons, which is usual for a monopole Hamiltonian 
[p^ . Deviations from this dependence are characteris¬ 
tic of the continuum coupling. For = 0 (a limit of 

one-neutron drip-line), one can see an odd-even stagger¬ 
ing (OES) of Ecorr around an average V-dependence. A 
blocking of the virtual scattering to the particle contin¬ 
uum by an odd nucleon diminishes the continuum cor¬ 
rection to the binding energy of broken-pair g.s. of odd- 
N nuclei. This ”drip-line effect”, which is seen also in 
the coordinate-space Hartree-Eock-Bogolyubov approach 
[p^ , is restricted to a narrow window of excitation ener¬ 
gies around = 0. It vanishes if Ecorr is calculated 

for = 4 MeV (see the dotted line with triangles in 

Eig. la). 

If is calculated in SMEC (the solid line), one 

finds an opposite OES effect with enhanced Ecorr for odd- 
N isotopes. This is due to which is lower in odd- 

N isotopes than in neighbouring even-V isotopes. This 
threshold induced OES effect leads to the attenuation 
of OES of binding energies. An enhanced modification 
of one-neutron separation energies Sn is seen for odd-V 
oxygen isotopes (see Eig. lb). In general, the effect of 
continuum coupling on one- and two-neutron separation 
energies is seen mainly for Sn <3 MeV, i.e. close to the 
one-neutron drip-line. 

In fluorine isotopes, both neutron-neutron (T = 1) and 
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neutron-proton (T = 0,1) couplings between states in 
Q and P are present. The nn-coupling has been ad¬ 
justed in the oxygen isotope chain and its value is kept 
unchanged. The optimal value of the strength of the np- 
continuum coupling : -|- /3), varies from 

a standard value ~ for well-bound isotopes 

at the valley of stability to near 

the neutron drip-line. This coupling modifies strongly 
T = 0 monopole terms of the SM Hamiltonian. For 
pd”p) ^ {1/2)Vq™\ good agreement for S'„, S 2 n and 
OES of binding energies : 


topes with those calculated in SMEC and in SM. SMEC 
results correspond to : = (1/2 )Eq^""'^ (the solid 

line) and (the dotted line). SM results 

have been obtained using two different SM interactions 
: (i) with the modified T = 0,1 monopole terms (the 
dashed line), as discussed above, and (ii) without these 
modifications (the dashed-dotted line). One can see that 
the SM calculations overestimate systematically the OES 
effect, in particular close to the neutron drip-line. The 
np-continuum coupling strongly attenuates the OES and 
may even wash it out for values of ratio 
close to the accepted value > 2 in nuclei from the valley 


= /2)[E{N+ l)-2E{N)+E{N-1)] (1) of stability. 


is found provided the T = 0 two-body interaction 
is modified as follows : = -1-0.6 MeV, 

<5En^=o =-h0.9MeV. 

005/2.0/7/2 

Fig. 2a shows the neutron number dependence of Ecorr 
for fluorine isotopes, which is calculated for Eg^"^^ = 
(l/2)kg^"'"^. Ecorr depends both on \Np — 7V„| and 
{Np + Nn) (compare Figs, la and 2a), where Np, Nn 
are the number of protons and neutrons in the valence 
space. Three cases shown in Fig. 2a can be directly com¬ 
pared with those shown in Fig. la. For En^’^'^ = 4 MeV, 
one can see the OES of Ecorr which is absent in the oxy¬ 
gen chain. This is the salient feature of the np-coupling. 
Ecorr in odd-odd isotopes is increased as compared to the 
neighboring odd-even ones. The size of Ecorr in this case 
is weakly dependent on the precise value of . Qual¬ 
itatively similar effect can be seen also at the neutron 
drip-line (see the results for = 0 in Fig. 2a) but 

now the OES due to the np-coupling is attenuated by 
an opposite staggering due to the nn-coupling (see also 
Eig. la). If En^’^^ is calculated in SMEC for each nu¬ 
cleus (the solid line), the OES of Ecorr around a smooth 
iV-dependence is enhanced due to combined effects of the 
np-continuum coupling, which weakly depends on , 
and the nn-continuum coupling, which closely follows the 
OES of These two effects act ’in phase’, hence 

contributing an additional binding for odd-V nuclei. 

A strong sensitivity of Ecorr in fluorine chain to the 
value of the strength can be seen in Eigs. 2b and 

2c which show the neutron number dependence of Sn and 
S 2 n- The dashed line shows the SM results, whereas the 
solid line gives SMEC results for Eg^”^^ = (l/2)Eg^””^ 
The agreement of SMEC results with the experimental 
data is excellent for neutron-rich nuclei. The dotted line 
in Figs. 2b and 2c corresponds to SMEC calculation 
with Eg*'"^^ = 2Eg*'""'^ Eor this standard value, one finds 
a somewhat better agreement with the experimental data 
for light fluorine isotopes. On the contrary, description 
of neutron-rich nuclei is poor. 

A useful indicator of OES is (Eq. |^). Fig. 3 

compares experimental data for A^^\N) in fluorine iso- 



FIG. 2. The same as in Fig. 1 but for the fluorine isotopes. 
SMEC calculations shown by the solid line are performed 
for the strengths of the continuum coupling : = 414 

MeV-fm^ and eJ"^' = (l/2)Eo^’""^ The dotted lines in (b) 
and (c) have been obtained for Eq"^^ = 2Eq""\ The closed 
circles with error bars shows the data ]l^. 
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In summary, we have performed first large-scale SM 
calculations including the coupling to the continuum 
states which extend from the valley of stability up to 
the neutron drip-line. Our analysis demonstrates an 
important renormalization of the effective SM Hamilto¬ 
nian by the continuum coupling; the principal effects of 
this renormalization depend on the position of the one- 
neutron emission threshold , on isospin T, its pro¬ 

jection Tz and on the number of valence particles. Cer¬ 
tain features of this renormalization can be simulated 
in SM studies by a suitable shifts of the two-body ma¬ 
trix elements. The dependence of iilcorr on both 
and the number of valence neutrons iV„ and protons Np 
implies, however, that the two-body monopoles in SM 
studies must contain the additional terms depending on 
En^'^\ \Np — Nn\, and {Np + Nn), which become impor¬ 
tant near the drip line. 



FIG. 3. The OES of binding energies (|^ for SM (the 
dashed line) and SMEC with different ratios of np- and 
nn- continuum coupling strengths : = (1/2 )Vq""^ 

(the solid line) and (the dotted line). The 

dashed-dotted line with open diamonds, shows SM results 
for the USD-I-KB interaction without the modification of 
T = 0,1 monopole terms discussed in the text. The closed 
points with error bars shows the data jl^ . 

We have found an important attenuation of the OES 
of binding energies close to the neutron drip-line by the 
np-continuum coupling. This attenuation does not mean 
that the neutron pairing correlations become weaker (see 
Fig. la) but that the np-continuum coupling ’falsifies’ 
AO)(jV) indicator of neutron pairing. 

The analysis of experimental data in SMEC indicates a 
signihcant attenuation of np-continuum coupling in nu¬ 
clei with large neutron excess. A strong sensitivity of 
Sn{N) and A^^\N) to the value of allows 

to estimate this attenuation quantitatively. 

The continuum coupling correction to energy of a SM 
state depends strongly on the position of this state with 
respect to En^’^^■ This dependence induce downward 


shifts of excitation energy for SM states close to En^’^\ 
The T^—dependence of these shifts (see Figs, la and 2a) 
implies that these SM states (e.g. intruder states) , even 
for a constant value of one-neutron emission threshold, 
will be pushed down by an amount which depends on Tz- 

The coupling to the two-neutron continuum has not 
been included in the present analysis. This coupling 
is expected to be less important if Sn < S 2 n, i-S. for 
IV < 18 in O- and F-isotope chains. Since is a 

smooth function of N and, moreover, varies slowly near 
the neutron drip-line, therefore we expect that the cou¬ 
pling to the two-neutron continuum produces a smaller 
effect on the OES of binding energies than the coupling 
to the one-neutron continuum. One should keep in mind, 
however, that additional attractive correlations between 
nucleons in the scattering continuum are responsible for 
the appearance of ’’Borromean systems” close to the one- 
neutron drip line [p^ . 

We thank W. Nazarewicz for valuable suggestions and 
F. Nowacki for useful discussions. 
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